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Abstract 

Modeling QCD at large temperature with a simple holographic five dimensional theory encoding 
minimal breaking of conformality, allows for the calculation of all the transport coefficients, up to 
second order, in terms of a single parameter. In particular, the shear and bulk relaxation times are 
provided. The result follows by deforming the AdS background with a scalar dual to a marginally 
relevant operator, at leading order in the deformation parameter. 
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1 Introduction and results 



The evolution of the Quark-Gluon Plasma (QGP) produced in heavy ion collisions is fairly 
well described, after a short period of thermalization and before hadronization, by relativistic 
hydrodynamics, consistently with the strongly coupled regime of such system [TJ. Numerical 
simulations of the hydrodynamic evolution of the QGP require as input the value of the 
transport coefficients. While recent simulations indicate that the evolution of the QGP 
should be quite insensitive to most of the second order coefficients (see for example [2j [3j H] ) , 
it is definitely sensitive to the value of the shear viscosity [5], and can be influenced in a 
sizable way by the bulk viscosity and possibly the relaxation times [31 HJ E]. Moreover, 
a complete characterization of the Quark-Gluon Plasma of QCD up to second order still 
requires the knowledge of the whole set of coefficients. 

There are currently no first-principle reliable calculations of almost all the second order 
coefficients for QCD at strong coupling: lattice results give some estimates of the viscosities 
and the shear relaxation time [7], but they are affected by considerable uncertainties (see 
for example [8]). In fact, actual simulations, lacking solid data for QCD, make often use, 
as benchmark values of the transport coefficients, of the ones derived from the gravitational 
dual of TV = 4 SYM [9] (in some cases together with the bound on the bulk viscosity proposed 
in [10] and a relation for the relaxation times from [IT]). While the M = 4 SYM values for 
the "shear" coefficients are expected to be in the right ballpark for QCD, they still concern 
an exactly conformal theory, and in particular the bulk viscosity and many of the second 
order coefficients are not determined. 

In order to improve this situation, the first step is to break conformal invariance. Since 
QCD is approximately conformal in the temperature window 1.5T C < T < 4T C , the confor- 
mality breaking effects can be treated perturbatively. In this situation, probably the simplest 
way of modeling QCD holographically is by a theory where conformality is slightly broken 
by a marginally relevant operator. 

The aim of this note is to point out that, in such a scenario, all the transport coefficients up 
to second order for the uncharged plasma are given in terms of a single parameter (weighting 
the conformality breaking) by making use of the results in pT2l fTT] : they are collected in Table 
[TJ In particular, the behavior of the shear and bulk relaxation times is briefly discussed in 
section II . 21 

There are surely more precise ways of modeling holographically QCD (none of which 
is of course completely correct). Nevertheless, the model considered in this note has the 
considerable advantage of full calculability, providing one of the few examples in which all 
the second order transport coefficients are determined. Moreover, the results in Table [TJ 
hold for any theory with gravity dual, where conformality is broken at leading order by a 
marginally (ir)relevant operator (dual to a scalar with the simplest possible potential (EJ)), 
including the cascading plasma [13] and the D3D7 plasmas [14J. 
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1.1 Notation 



Uncharged relativistic hydrodynamics is determined, up to second order in the derivative 
expansion, by seventeen transport coefficients, fifteen of which are possibly independent [9], 
[TTj . On a general space with metric g^, the energy momentum tensor 

T*» = eu^vT + pA^ + tt^ + A^n , where A^-f + (1) 

is determined by the energy density e, fluid velocity u M {u^u^ = — 1), the transport coeffi- 
cients in its "viscous shear" part: 



71-"" = -rja^ + rjr^ 



V • 7/ i r 

(DO + o^ v + k R <lxl/> - 2u a u p R a< » v>p 
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+A 2 a< / ^ >A + A 3 fi <M A ^T >A + K*2u a upR a<,xu>f) 

+r]rl^—a txv + A 4 V<" log sV^ log s (2) 
and in its "viscous bulk" part: 

n = -C(V- M ) + CrnD(V- M )+6^^ + 6(V-n) 2 + ^^ + ^V^log S V^log S 

+^R + ^u a u"R a(5 , (3) 

while the pressure is given by the equation of state p(e). The various structures in (j2J) and 
([3]), apart from the obvious Riemann and Ricci tensors and scalar curvature (R^ upa , R^ u ', R), 
are given by: 
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D = m m V m , V^A^V,, (J" 1 ' = V>" + - -A^(V • u) 



^ = ~(vy- v^o, (4) 

and for a generic tensor it was used the notation: 

(A"") = A<^> = l -A^A^{A aP + A Pa ) - l -A^A^A af} . (5) 

Finally, s is the entropy density, while the speed of sound will be denoted as cj. = dp/de. 

The shear viscosity r\ and the second order coefficients T n ("shear" relaxation time), k, 
Ai, A2, A3 are the only ones defined in conformal fluids, as the one of M = 4 SYM. All the 
others coefficients, i.e. the bulk viscosity ( and the second order coefficients k*,t*,\^,tu 
("bulk" relaxation time), £1, £2, £3, £4, £5, £e> are only defined in non-conformal plasmas. 



1.2 The estimate 

Consider a gravity dual model for QCD at large temperature, where the leading conformality 
breaking effect is captured by adding to the five dimensional metric a non trivial dilaton 
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profile, dual to a marginally relevant operator. The main observation of this note is that, for 
the simplest scalar potential, the transport coefficients are completely determined in terms 
of a single parameter. Defining: 

5 = (1 - 3c 2 s ) , (6) 
at first order in 5 the transport coefficients are given in Table [U This result, which is the 
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Table 1: The transport coefficients, in the notation of ([I])-®, for a marginally (ir)relevant 
deformation of a conformal theory, at leading order in the deformation parameter 8 = (1 — 
3c^). The holographic equation of state is e = 3(1 + 5)p. 

main content of this note, follows directly from [12], [TT] (which already contains a part of 
the relations in Table [Jl]) and will be derived in section [2j 

Possibly the main novel results contained in Table [1] concern the two relaxation times 
r n , tjj. Specifically, at leading order in the conformality breaking, the bulk relaxation time 
Tn is not proportional to the bulk viscosity. The behavior of the shear relaxation time is 
instead more interesting, since it depends on the speed of sound. For a phenomenologically 
realistic behavior of the latter, t w is decidedly increasing when reducing the temperature. In 
particular, it increases faster than r n . 

^ee also [in Hoi ng mi Ha ng na ^ . 
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Moreover, using the above results it is easy to verify that the relation 

4Ai + A 2 = 2777V, (7) 

holds, at first order in 5. It has been shown in [2U E2] that ((?]) is satisfied in all the known 
examples of conformal plasmas (in d > 4 spacetime dimensions, with of without conserved 
global charges) with dual gravity description. Our results provide a unique validity check of 
the above relation in non-conformal settingsj§ 

While the present system can model at best the regime of QCD away from the critical 
temperature, where there are certainly other ways of modeling QCD, it would be unexpected 
if the qualitative behavior of the transport coefficients described above turned out to be 
drastically different. 

In order to give an illustrative example of numerical estimates of the coefficients, we have 
to chose one input parameter. As in [4], we use the results for the speed of sound from the 
lattice study in [23J. We consider a temperature T ~ 1.5T C , which is a reasonable value for 
the RHIC experiment. Then from [23] we read c 2 s ~ 0.283 from which we get the numbers 
in Table EJ The reported values provide corrections up to 17% to the conformal ones (when 
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Table 2: The transport coefficients at T ~ 1.5T C and c 2 s ~ 0.283. 
2 We thank Todd Springer for this observation. 
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the latter are defined). In particular 2ttTtk = 1.394 is a bit larger than the conformal value 
(1.307) and more similar to the one used in [3]. The numerical difference is by definition not 
very large, but sizable. Obviously, increasing the temperature reduces this difference and at 
T ~ 3T C , which could be a significant temperature for LHC, the corrections to the conformal 
values are below 10%. 

Let us conclude this section by describing the approximations involved in applying these 
relations to QCD. First of all, QCD does not have a purely gravitational dual. Nevertheless, 
experience teaches that simple gravity models of metric plus scalar are in good quantitative 
agreement with the lattice results for certain observables. In particular, we are interested 
in the regime, relevant in the early stages of the QGP evolution, at T > T c where QCD is 
nearly conformal and strongly coupled. Moreover, in the hydrodynamic regime the gravity 
description and actual QCD are in good agreement (e.g. the result for the shear viscosity). 
On the other hand, in QCD the gluon condensate is marginally relevant in the asymptotically 
free regime, while in the experimental regime we are interested in, the theory is strongly 
coupled and this operator can be expected to have developed a sizable anomalous dimensions 
The other caveat concerns the effects of the flavors and the chemical potential, which are 
not accounted for in Table HJ but are expected to give the latter subleading corrections. 

In view of these considerations, the relations in Table [1] can provide a fair estimate^ of the 
initial behavior of the hydrodynamic evolution at RHIC and LHC. 

2 Derivation 

The leading conformality breaking effects of a source for a marginally (ir)relevant operator 
can be captured in the dual gravitational setting by a so-called Chamblin-Reall model. 
Consider an effective five dimensional theory with metric plus a single scalar with potential 
V (</>)• It models the breaking of conformality at leading order in a small parameter e if 
V(<f))\ e= o = Vo, where the negative cosmological constant V allows for and AdS solution. 
The operator dual to (on the unperturbed AdS solution) is of dimension four, that is it is 
marginally (ir)relevant, if 9|y(0)|^ =o = 0(e 1+a ) with positive a. The simplest such class of 
models, and the one we are interested in, is given by: 



3 This situation could be modeled with a scalar dual to a relevant operator [TBI US]- I n this case the 
computation of the second order coefficients is highly more complicated. 
4 Better than the one provided by TV = 4 SYM [9]. 



V{cj ) )=Vo + 6 ( t ) + 0(6 1+a ). 




At leading order: 



V{<j>) ~ V e^ v ° , 



(9) 
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i.e. the model is in the Chamblin-Reall class [21 

For this class of models, the proof of the relations in Table [1] follows directly from [12] . 
Let us summarize it. The starting point is the fact that Chamblin-Reall models in d + 1 
dimensions, for particular values of the coefficient of the exponential in the potential, can 
be obtained from dimensional reduction on a 2a — d torus of pure gravity plus cosmological 
constant in 2a + 1 dimensions. This happens when the parameter a, which determines 
together with d the coefficient in the exponential, is semi-integer. For these values of a, one 
can then start from the well-known AdS^a+i solution and its dual hydrodynamic energy- 
momentum tensor, and obtain the hydrodynamic energy-momentum tensor for the dual to 
the Chamblin-Reall model by simple toroidal dimensional reduction. 

The crucial observation in [12] is that, from the point of view of the theory in d + 1 
dimensions, the equations are smooth in the parameter a. This allows for the computation 
of the hydrodynamic energy-momentum tensor for arbitrary values of a > d/2% 

The procedure is as follows. One starts from a Chamblin-Reall model in d + 1 dimensions 
for whatever a > d/2 and performs the continuation (which is smooth) to the nearest value 
a which is semi-integer. The latter theory is the compactification of a theory admitting 
a AdS2d+\ solution, so its dual energy-momentum tensor, which will be a function of a, 
can be calculated straightforwardly. This energy-momentum tensor can thus be continued 
(smoothly) back to the one of the theory corresponding to the original value a. 

In particular, all the transport coefficients for this theory will automatically be determined 
by the conformal ones in the higher dimensional theory, modulo an overall constant (the 
volume of the torus) which can be fixed knowing just one coefficient. 

Let us see concretely how this procedure is implemented. One can determine a by the 
fact that the equation of state in these models is e = (2cr — l)p [12], so that a = 2 + 35/2 
in the notation of Table [TlFI Thus, for a small deformation S of a conformal theory, a = 5/2 
and the relevant starting solution is AdS2a+i = AdS^, whose dual conformal hydrodynamics 
was considered in [25]. Let us write the results of [25J in the present notation: 
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5 To be precise, with unit AdS radius, Vb = — 12 and, from the calculation of the speed of sound, e 2 = 965 

6 At a = d/2 the action is singular [12] . 
7 And cr = 2 + e 2 /64 in the notation of ©. 
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The procedure to obtain the desired coefficients involves reducing the energy momentum 
tensor on a circle (2a — d = 1) of volume V, continuing it back to a and expanding it at 
first order in S p2] ; examples of results of this procedure are (the arrows denote the analytic 
continuation): 
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1 El- The other coefficients 



(2a-l)(d-l) 

where the leading "conformal" term in ( I13p fixes the value V 
in Table [1] are obtained in the same way. 

Let us conclude by stressing again the fact that the relations in Table [T] are valid for any 
theory where conformality is broken at leading order by a marginally (ir) relevant deforma- 
tion, with the dual scalar having the potential ©. These theorieqfl include the cascading 
plasmas [T3] and the D3D7 plasmas [J4]. The latter are the first examples of holographic 
plasmas including the effects of dynamical flavors in a completely controllable framework. 
In this case, the relations in Table [1] match precisely the coefficients calculated in |2oJ^| and 
complete the determination of all the second order transport coefficients in those plasmas. 
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